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X-ray Observations of Radio Transients without Optical Hosts 

Steve Croft 1 , John A. Tomsick 2 , and Geoffrey C. Bower 1 

ABSTRACT 



We present a 50 ks Chandra ACIS-I X-ray observation of the Bower et al. VLA 
£vq archival field. The observations reach a limiting sensitivity of ~ 10 -4 counts s _1 , 

corresponding to a flux of a few times 10 _15 erg s _1 cm -2 for the models we 
h- 5 explore. The Chandra observations were undertaken to search for X-ray coun- 

terparts to the eight transient sources without optical counterparts, and the two 



transient sources with optical counterparts seen by Bower et al. Neither of the 
sources with optical counterparts was detected in X-rays. One of the eight op- 
tical non-detections is associated with a marginal (2.4cr) X-ray detection in our 



Chandra image. A second optically-undetected |Bower et al.| transient may be 



associated with a z = 1.29 X-ray detected quasar or its host galaxy, or alterna- 
tively is undetected in X-rays and is a chance association with the nearby X-ray 

source. The X-ray flux upper limits, and the one marginal detection, are con- 
ed 



sistent with the interpretation of Ofek et al. that the optically-undetected radio 



transients are flares from isolated old Galactic neutron stars. The marginal X-ray 
>■ detection has a hardness ratio which implies a temperature too high for a simple 

one-temperature neutron star model, but plausible multi-component fits are not 
excluded, and in any case the marginal X-ray detection may be due to cosmic 
rays or particle background. The X-ray flux upper limits are also consistent with 
Q flare star progenitors at > 1 kpc (which would require the radio luminosity of 

the transient to be unusually high for such an object) or less extreme flares from 
brown dwarfs at distances of around 100 pc. 
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Introduction 



In recent years, authors such as Bower et al. (2007), hereafter B07; Becker et al. (2010); 



Bower & Saul ( |2011 ); Bannister et al. (2011); Bell et al. (2011); and others, have made use of 
archival data (observations of calibration fields, or in some cases, entire telescope archives) 
to search for radio transients. Other studies have targeted particular fields with observations 



dedicated to the search for transients (Ofek et al. 2011, and references therein). 



Such studies are preparing the groundwork for a new generation of high-throughput 
radio telescopes: both dedicated new observatories - - for example, the Allen Telescope 



Array (Welch et al. 2009), LOFAR (Rottgering 2003), and ASKAP (Johnston et al. 2008) 



and upgrades to older instruments -- for example, Apertif (Verheijen et al. 2008) and 



EVLA (Perley et al. 2009). A major goal for these new facilities is to survey large areas of 



sky at high sensitivity, both to build up deep, wide-field images, and to search for transient 



and variable sources from a range of progenitors (Lazio et al. 2008). 



Currently, large areas of rate versus sensitivity parameter space remain to be explored 



(Bower et al. 2010 Croft et al. 2011), and the radio transient population is not well un- 



derstood. In observations where radio transients are detected, sometimes counterparts are 
seen in existing surveys at other wavelengths. For example, transients may coincide with the 
nuclear regions of galaxies, in which case (depending on the galaxy redshift and hence the 
inferred radio luminosity of the transients) they are most likely due to active galactic nuclei 
(AGNs), or perhaps supernovae (SNe) or gamma ray bursts (GRBs). Or they may be offset 
from the peak of the optical emission, as is the case for two of the transients seen by B07, 
leading to their interpretation as SNe or GRBs. 

For some radio transients, however, no counterparts are seen in archival images at other 
wavelengths, and no follow-up observations taken close in time to when the transients were 
seen are available. Follow-up observations will become more common as radio observatories 
move increasingly into a regime of issuing transient alerts for rapid follow-up by other facili- 
ties, but the field is still in its infancy compared to optical transient facilities which are now 



issuing transient alerts as a matter of routine (e.g., Law et al. 2009). However, the lack of 



galaxy counterparts in deep optical images can help exclude some models for these sources 
(such as AGNs, SNe and GRBs, unless these sources are at very high redshift or extremely 
obscured), and the lack of stellar counterparts can also constrain models with Galactic pro- 
genitors for these events. B07 consider radio SNe, GRBs, late-type stars, soft gamma-ray 
repeaters, X-ray binaries, pulsars, microlensing, and reflected solar flares as possible progen- 
itors. They conclude that late-type stars are the most likely progenitors for transients with 
mJy radio fluxes and no optical counterparts. 



Ofek et al. (2010), hereafter O10, suggest old isolated Galactic neutron stars as a pos- 



sible progenitor for the sources without optical identifications, and Nakar & Piran (2011) 
suggest that the sources with optical counterparts may be due to neutron star - neutron star 
or neutron star - black hole mergers. The latter are of particular interest because they are 
expected to also be producers of gravitational waves. X-ray observations can be a powerful 
discriminant among progenitor models, but information on the X-ray properties of the tran- 
sients and transient hosts is limited. With this in mind, we obtained Chandra observations 
of the B07 field. 



We assume an Vt m = 0.27, Vt\ = 0.73, H = 71 km s 1 Mpc 1 cosmology (Jarosik et 



al. 2011). Magnitudes are given in the Vega system unless stated otherwise. 



Observations 



2.1. Radio, optical, and Infra-Red Observations of the Bower et al. (2007) field 



B07 analyzed archival data at 5 and 8.4 GHz of a Very Large Array (VLA) calibration 
field centered at 15 h 02 m 20!53, 78° 16' 14V905. The data spanned ~ 20 yr with 944 epochs, 
approximately one per week, each with an integration time of ~ 20 min, all on the same field. 
The images had a full-width at half-maximum (FWHM) of 8'6 and 5'1 at 5 and 8.4 GHz 
respectively, and typical root-mean-square (rms) fluctuations at the center of the images 
were ~ 50 /iJy per epoch. 

B07 compared the images from epoch to epoch, and detected 7 single-epoch transient 
sources at 5 GHz, and 1 at 8.4 GHz. The transients are listed in Table [1} B07 denoted 
these sources with identifiers beginning "RT" and the epoch in which they were observed; 
for convenience and clarity here we denote the 5 GHz and 8 GHz transients with identifiers 
beginning "5T" and "8T" respectively, although their formal identifiers are those beginning 
"RT" with an epoch of observation as cross-referenced in Table [Tj For each of these transients, 
there was no simultaneous coverage at the other of the two observing frequencies, and no 
detection of the sources in longer-timescale average images. 
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B07 also constructed a deep image of the field using all data at each frequency. The 
resulting images at 5 and 8.4 GHz had an rms of 2.6 and 2.8 /xJy, respectively. Sources which 
were seen in many of the individual epochs and also in the deep field image were termed 
"steady" by B07. They saw 8 steady sources at 5 GHz, shown in Table [I] with identifiers 
beginning "5S", and 4 of these correspond to steady sources in the deep image at 8.4 GHz, 
shown with labels beginning "8S" in Table [TJ 

B07 also compared images made by making images from two months' worth of data at 
a time, and found two additional "long-timescale" transients, one at each frequency, shown 
with identifiers beginning "5L" and "8L" in Table [TJ 

Optical and infra-red (OIR) observations of the field were also reported by B07, both 
with LRIS on Keck I, to limiting magnitudes of g ~ 27.6 mag and R ~ 26.6 mag, and with 
PAIRITEL to J « 19.2 mag, H fa 18.5 mag and K s fa 18.0 mag. Of the two long-duration 
transients, one has a match with a, z = 0.25 galaxy. Of the eight short-duration transients, 
one has a likely match (with a. z = 0.04 galaxy). O10 also present observations in K s -hand 
using the Hale 5.08-m telescope to a limiting magnitude of K s ranging from 19.2 to 20.4 
for those sources undetected in OIR by B07. They found no additional identifications. Two 
transients are referred to by B07 as "possible" matches, one with a galaxy at z — 0.25, and 
the other with a galaxy with unknown redshift (see Table 111), but the images of O10 suggest 
that both of these associations are unlikely. We consider these two, like the remaining six 
B07 radio transients, to be unidentified in all the OIR images from B07 and O10. 

As discussed above, B07 suggest a range of possible progenitors for these transients, 
and O10 expand on this discussion. O10 also suggest that the transients seen at the Nasu 



Pulsar Observatory by Matsumura et al. (2009) and other papers from the same group may 



be members of the same class. Recently, Croft et al. (2010), Croft et al. (2011), and Bower 



& Saul (2011) suggested that the Nasu transients are probably not astronomical events, 



and that their quoted rate is at odds with other measurements, although the parameters of 



the Nasu survey are not clear from the literature (see Ofek et al. 2011). However, the B07 



transient rate is consistent with measurements from other surveys (see, for example, Croft 



et al. 2011, Fig. 12), and two of the B07 transients have clear OIR counterparts. O10 show 



that for most classes of object, in particular extragalactic events including supernovae and 
gamma ray bursts, one would expect a large fraction of the transients to have associated 
host galaxies. They argue that the most attractive explanation is that the transients are 
associated with isolated old Galactic neutron stars. In order to test this hypothesis, and the 
flare star scenario discussed by B07, we obtained X-ray observations of the field. 



2.2. X-Ray Observations 

X-ray observations were obtained with Chandra during Cycle 11, using the ACIS-I 
instrument with no grating or filter. The exposure time was 49,441 s and the observation 
start date was UT 2010 July 3 01:14:32. The observation was centered at 15 h 01 m 50:58, 
78° 16' 50 / . / 90 so as to encompass the majority of the VLA field. 

For our analysis, we started with the "level 2" data products, and these were processed 
at the Chandra X-ray Center with pipeline ("ASCDS") version 8.3. We performed fur- 
ther processing with the Chandra Interactive Analysis of Observations (CIAO) version 4.3 
software and version 4.4.1 of the Calibration Data Base (CALDB). For the work described 
below, we used an image from the four ACIS-I CCDs, including photons in the 0.3 — lOkeV 
bandpass. 

A catalog was created from the full-band image, by running the CIAO source detection 
algorithm "wavdetect". We used 4 different pixel binnings (1, 2, 4, and 8) in order to detect 
sources that are either truly extended or sources that appear extended because the PSF 
broadens off-axis. We combined the catalogs from these 4 runs. A total of 128 sources were 
detected, and their properties are shown in Table [2j In each run, the threshold was set such 
that 1 spurious source would be expected, so one would expect up to 4 of the sources in 
Table [2] to be spurious. In most cases, the parameters shown in Table [2] correspond to the 
lowest pixel binning in which a source was detected. In some cases (such as the extended 
source X29, or off-axis sources), the lowest binning clearly did not enclose all of the counts 
and a higher binning was used. The binning corresponding to the reported source parameters 
is shown in Table [2j 



Table 2. Chandra sources 



Chandra ID 


Chandra RA 


Chandra Dec 


Chandra pos. 


Bin b 


ACIS counts c 


Hardness VLA 




(J2000) 


(J2000) 


unc. (arcsec) a 




(0.3- lOkeV) 


constraint" 1 ID C 


XI 


14 57 57.57 


78 15 29.7 


1.46 


1 


45.4 ±9.5 


< 0.71 


X2 


14 58 14.94 


78 16 38.1 


1.63 


1 


60.0 ± 12.2 


< 0.89 


X3 


14 58 31.24 


78 18 14.5 


1.25 


2 


28.6 ±6.8 


< 0.04 


XI 


14 58 33.49 


78 15 03.8 


1.29 


2 


45.4 ±8.9 


-0.29 ±0.22 


X5 


14 58 45.15 


78 15 29.3 


0.95 


1 


8.0 ±3.2 


< 0.62 


X6 


14 58 46.87 


78 18 60.0 


1.98 


1 


23.9 ±7.8 




X7 


14 58 48.87 


78 15 13.5 


0.86 


1 


7.9 ±3.0 


< 0.26 


X8 


14 59 02.72 


78 16 58.6 


1.22 


3 


22.8 ±6.4 


0.03 ± 0.40 


X9 


14 59 11.82 


78 20 50.5 


1.11 


3 


36.1 ±8.1 


-0.33 ±0.32 


X10 


14 59 20.11 


78 17 05.8 


0.85 


1 


8.1 ±3.2 


< 0.07 


Xll 


14 59 23.81 


78 18 08.5 


1.02 


2 


9.8 ±3.7 


< 0.31 


X12 


14 59 25.15 


78 15 20.7 


1.40 


3 


30.4 ±8.0 


0.24 ±0.41 


X13 


14 59 26.92 


78 14 28.5 


1.37 


3 


38.6 ±8.0 


0.36 ±0.29 


X14 


14 59 26.98 


78 16 49.8 


0.93 


1 


25.2 ±6.1 


0.09 ±0.25 


X15 


14 59 30.92 


78 14 05.9 


1.28 


3 


22.8 ±6.6 


< 0.22 


X16 


14 59 32.13 


78 11 20.2 


1.69 


1 


24.8 ±8.1 




X17 


14 59 35.61 


78 12 59.4 


0.93 


2 


23.7 ±5.8 


-0.60 ±0.31 


X18 


14 59 45.84 


78 21 09.6 


1.01 


2 


16.6 ±4.8 


0.32 ±0.43 


X19 


14 59 46.99 


78 20 23.7 


0.83 


2 


40.5 ± 7.2 


-0.20 ±0.20 5T7 


X20 


14 59 50.62 


78 10 46.7 


1.43 


3 


12.9 ±5.0 


0.14 ±0.75 


X21 


14 59 55.38 


78 11 16.7 


1.16 


3 


24.8 ±6.9 


0.30 ±0.41 


X22 


14 59 57.59 


78 13 49.5 


0.72 


1 


9.6 ±3.5 


-0.28 ±0.51 


X23 


14 59 57.89 


78 21 27.6 


0.83 


1 


26.8 ±6.2 


-0.36 ±0.25 


X24 


15 00 00.57 


78 22 53.2 


0.81 


1 


8.4 ±3.3 


< 0.31 


X25 


15 00 02.85 


78 17 41.3 


1.02 


2 


15.5 ±4.6 


> -0.03 


X26 


15 00 05.30 


78 14 47.3 


0.70 


1 


7.5 ±2.8 


< 0.21 


X27 


15 00 06.04 


78 18 16.1 


0.74 


1 


6.6 ±2.8 


-0.21 ±0.54 


X28 


15 00 06.45 


78 18 35.1 


0.69 


1 


48.2 ±7.7 


-0.07 ±0.16 


X29 


15 00 08.62 


78 12 48.6 


1.16 


1 


166.8 ±21.1 


-0.62 ±0.11 5S1 


X30 


15 00 09.28 


78 16 18.9 


0.82 


1 


18.1 ±4.7 


0.40 ± 0.33 


X31 


15 00 12.56 


78 21 03.4 


1.11 


1 


16.1 ±5.0 


-0.45 ±0.46 


X32 


15 00 13.00 


78 15 45.8 


1.04 


2 


16.4 ±4.9 


> 0.24 


X33 


15 00 15.41 


78 24 20.9 


1.38 


2 


18.7 ±5.9 


-0.38 ±0.59 


X34 


15 00 19.57 


78 14 04.8 


1.00 


2 


12.5 ±4.1 


-0.49 ±0.48 


X35 


15 00 24.61 


78 13 41.5 


0.71 


1 


30.0 ±5.8 


-0.30 ±0.22 


X36 


15 00 25.27 


78 15 49.5 


0.77 


1 


8.5 ±3.2 


> -0.07 


X37 


15 00 26.05 


78 15 18.1 


0.73 


1 


23.5 ±5.1 


-0.39 ±0.31 


X38 


15 00 27.14 


78 11 29.0 


1.43 


3 


15.1 ±5.4 


-0.02 ±0.62 


X39 


15 00 29.41 


78 13 02.9 


0.76 


1 


12.6 ±4.0 


-0.31 ±0.31 


X40 


15 00 32.58 


78 21 12.0 


0.90 


2 


9.6 ±3.6 


> -0.54 


X41 


15 00 33.99 


78 14 58.5 


0.76 


1 


6.3 ±2.6 


< 0.59 


X42 


15 00 35.55 


78 21 36.6 


0.66 


1 


332.7 ± 19.9 


-0.37 ±0.06 


X43 


15 00 43.05 


78 09 38.2 


1.06 


2 


16.5 ±4.9 


-0.48 ±0.47 


X44 


15 00 44.19 


78 26 41.8 


1.32 


1 


103.1 ± 14.3 


0.05 ±0.19 


X45 


15 00 45.54 


78 12 12.1 


0.71 


1 


4.7 ±2.2 





Table 2 — Continued 



Chandra ID 


Chandra RA 


Chandra Dec 


Chandra pos. 


Bin b 


ACIS counts c 


Hardness VLA 




(J2000) 


(J2000) 


unc. (arcsec) a 




(0.3-10keV) 


constraint d ID C 


X46 


15 00 47.01 


78 21 20.2 


0.86 


2 


43.5 ±7.7 


0.09 ±0.18 


X47 


15 00 48.72 


78 18 00.3 


1.12 


3 


19.4 ±6.2 


> 0.16 


X48 


15 00 52.64 


78 09 00.6 


0.85 


3 


65.2 ±9.3 


-0.38 ±0.18 


X49 


15 00 54.07 


78 14 17.9 


0.97 


2 


3.7 ±2.2 




X50 


15 00 54.39 


78 08 13.2 


1.31 


3 


14.9 ±5.2 


< 0.02 


X51 


15 00 55.19 


78 13 48.3 


0.72 


1 


26.2 ±5.6 


-0.45 ±0.27 


X52 


15 00 55.43 


78 17 56.8 


0.67 


1 


4.5 ±2.2 


< 0.26 


X53 


15 00 58.19 


78 15 11.7 


1.55 


4 


36.0 ±10.6 


0.23 ±0.49 


X54 


15 01 05.65 


78 12 56.3 


0.72 


1 


21.0 ±4.8 


0.13 ±0.34 


X55 


15 01 05.82 


78 15 39.6 


0.76 


1 


6.0 ±2.6 


< -0.06 


X56 


15 01 15.61 


78 23 02.7 


0.81 


1 


35.8 ±7.2 


-0.07 ±0.19 


X57 


15 01 20.47 


78 20 29.1 


0.80 


1 


7.6 ±3.0 


< 0.05 


X58 


15 01 22.67 


78 18 06.1 


0.68 


1 


38.5 ±6.4 


0.52 ± 0.23 5S5/8S2 


X59 


15 01 23.82 


78 09 55.6 


0.91 


2 


16.9 ±4.8 


< -0.23 


X60 


15 01 23.88 


78 06 47.9 


0.73 


2 


165.5 ± 15.1 


-0.39 ±0.09 


X61 


15 01 26.36 


78 24 39.6 


0.90 


1 


33.3 ±7.5 


-0.16 ±0.20 


X62 


15 01 26.86 


78 07 38.9 


1.37 


3 


14.7 ±5.5 


0.09 ± 0.45 


X63 


15 01 29.75 


78 21 19.3 


0.72 


1 


17.0 ±4.5 


-0.08 ±0.40 


X64 


15 01 29.75 


78 19 14.0 


0.73 


1 


7.0 ±2.8 


0.27 ±0.58 


X65 


15 01 30.50 


78 25 11.9 


0.78 


2 


109.2 ± 13.2 


-0.41 ±0.13 


X66 


15 01 31.14 


78 20 18.4 


0.75 


1 


14.1 ±4.0 


0.34 ± 0.55 


X67 


15 01 32.91 


78 12 46.6 


1.00 


2 


13.0 ±4.2 


< -0.10 


X68 


15 01 35.97 


78 06 53.2 


1.04 


3 


32.2 ±7.7 


-0.32 ±0.32 


X69 


15 01 37.58 


78 10 03.2 


0.76 


1 


6.2 ±2.6 




X70 


15 01 41.14 


78 09 38.5 


1.28 


3 


12.5 ±4.8 


< 0.74 


X71 


15 01 44.79 


78 15 26.6 


0.64 


1 


172.6 ± 13.3 


-0.25 ±0.09 


X72 


15 01 45.39 


78 08 36.8 


0.77 


1 


80.1 ±10.5 


-0.61 ±0.14 


X73 


15 01 45.95 


78 12 03.9 


0.73 


1 


24.3 ±5.2 


0.53 ±0.33 


X74 


15 01 48.39 


78 15 41.0 


1.04 


2 


7.0 ±3.0 




X75 


15 01 53.42 


78 05 34.5 


1.64 


4 


18.4 ±5.7 


> -0.27 


X76 


15 01 54.60 


78 09 08.5 


1.06 


1 


11.4 ±3.9 


< 0.37 


X77 


15 01 54.69 


78 26 13.5 


1.83 


4 


22.1 ±7.6 


< 0.59 


X78 


15 01 56.25 


78 18 59.5 


1.72 


4 


24.0 ±9.0 


< 0.48 


X79 


15 01 59.30 


78 10 35.0 


0.68 


1 


6.7 ±2.6 




X80 


15 02 02.31 


78 15 22.1 


0.66 


1 


41.9 ±6.7 


-0.01 ±0.18 


X81 


15 02 04.83 


78 08 47.5 


0.78 


2 


5.2 ±2.6 




X82 


15 02 06.85 


78 05 16.9 


2.28 


4 


33.1 ±8.4 


> 0.04 


X83 


15 02 07.00 


78 16 22.7 


0.66 


1 


46.2 ±6.9 


-0.32 ±0.21 


X84 


15 02 07.50 


78 15 48.7 


0.68 


1 


19.7 ±4.6 


-0.33 ±0.33 


X85 


15 02 08.80 


78 16 32.3 


0.72 


1 


8.7 ±3.2 


-0.11 ±0.49 


X86 


15 02 12.95 


78 09 08.6 


1.37 


3 


13.1 ±4.7 




X87 


15 02 14.95 


78 18 30.6 


0.68 


1 


5.4 ±2.4 


< 0.59 


X88 


15 02 15.93 


78 11 38.1 


0.73 


1 


11.2 ±3.6 


-0.26 ±0.39 


X89 


15 02 17.58 


78 08 36.9 


0.89 


1 


9.9 ±3.6 


< -0.26 


X90 


15 02 19.21 


78 14 47.7 


0.72 


1 


14.6 ±4.0 


-0.02 ±0.41 
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Table 2 — Continued 



Chandra ID 


Chandra RA 


Chandra Dec 


Chandra pos. 


Bin b 


ACIS counts c 


Hardness VLA 




(J2000) 


(J2000) 


unc. (arcsec) a 




(0.3- lOkeV) 


constraint" 1 ID C 


X91 


15 02 21.11 


78 24 34.5 


1.53 


4 


44.1 ±11.3 


> 0.08 


X92 


15 02 24.68 


78 16 52.0 


0.72 


1 


14.5 ±4.0 


0.11 ±0.42 


X93 


15 02 24.92 


78 13 48.9 


1.34 


2 


13.5 ±4.4 


> 0.16 


X94 


15 02 27.06 


78 08 42.9 


0.92 


2 


8.5 ±3.3 


< 0.40 


X95 


15 02 35.73 


78 23 48.0 


0.89 


1 


32.7 ±7.0 


-0.38 ±0.22 


X96 


15 02 35.97 


78 10 51.2 


1.04 


2 


11.7 ±4.1 


-0.26 ±0.53 


X97 


15 02 37.60 


78 16 27.0 


0.67 


1 


3.7 ±2.0 




X98 


15 02 38.98 


78 09 02.3 


0.96 


1 


20.5 ±5.7 


-0.18 ±0.30 


X99 


15 02 41.57 


78 23 43.5 


1.02 


3 


24.8 ±6.4 


-0.35 ±0.34 


X100 


15 02 44.86 


78 20 04.3 


0.68 


1 


5.4 ±2.4 


< 0.26 


X101 


15 02 46.47 


78 17 56.7 


0.69 


1 


11.9 ±3.6 


-0.03 ±0.53 


X102 


15 02 47.44 


78 12 36.0 


1.02 


2 


22.9 ±5.8 


0.25 ±0.33 


X103 


15 02 48.22 


78 07 59.0 


1.58 


1 


24.2 ±8.2 


> -0.11 


X104 


15 02 51.38 


78 17 35.6 


0.70 


1 


32.1 ±5.9 


-0.01 ±0.23 


X105 


15 02 51.80 


78 11 27.0 


0.81 


1 


10.9 ±3.6 


-0.04 ±0.62 


X106 


15 02 53.37 


78 09 30.5 


0.96 


2 


19.3 ±5.3 


0.21 ±0.39 


X107 


15 02 56.01 


78 10 21.6 


1.25 


3 


26.7 ±6.9 


0.32 ±0.36 


X108 


15 03 02.99 


78 18 20.6 


0.70 


1 


26.2 ±5.4 


-0.12 ±0.27 


X109 


15 03 03.61 


78 14 51.9 


0.76 


1 


9.2 ±3.3 


-0.07 ±0.36 


XllO 


15 03 04.48 


78 22 13.0 


1.55 


4 


21.9 ±7.5 


< 0.41 


Xlll 


15 03 07.02 


78 17 48.9 


0.71 


1 


20.7 ±4.9 


-0.22 ±0.26 


X112 


15 03 10.12 


78 21 36.5 


1.88 


4 


31.2 ±8.9 


> -0.45 


X113 


15 03 11.45 


78 18 48.5 


0.74 


1 


9.7 ±3.5 


< -0.20 


X114 


15 03 13.38 


78 17 47.0 


1.22 


2 


13.8 ±4.4 


0.28 ± 0.43 


X115 


15 03 14.18 


78 22 51.3 


1.44 


3 


39.1 ±7.9 


< 0.43 


X116 


15 03 16.02 


78 10 04.1 


1.86 


4 


29.0 ±9.6 


> 0.53 


X117 


15 03 16.31 


78 13 02.8 


0.78 


1 


17.9 ±4.6 


-0.43 ±0.34 


X118 


15 03 23.14 


78 14 28.4 


0.68 


1 


97.6 ± 10.6 


-0.22 ±0.11 


X119 


15 03 24.06 


78 18 08.4 


1.34 


4 


21.7 ±8.9 


0.38 ±0.52 


X120 


15 03 25.76 


78 14 39.9 


0.75 


1 


13.8 ±4.1 


-0.29 ±0.43 


X121 


15 03 29.45 


78 16 06.9 


0.81 


1 


16.3 ±4.4 


< -0.20 


X122 


15 03 37.93 


78 18 14.9 


0.95 


2 


11.4 ±3.9 


0.38 ± 0.44 


X123 


15 03 42.10 


78 17 08.3 


0.72 


1 


15.2 ±4.2 


0.07 ±0.29 


X124 


15 03 45.73 


78 15 40.7 


0.95 


2 


8.8 ±3.6 


0.38 ± 0.44 8L1 


X125 


15 04 07.84 


78 13 16.4 


1.50 


4 


39.9 ± 10.3 


> -0.11 


X126 


15 04 40.84 


78 15 33.8 


1.96 


4 


40.2 ±10.4 




X127 


15 04 56.27 


78 16 42.8 


1.53 


4 


37.9 ±9.8 


> 0.31 


X128 


15 05 13.24 


78 18 17.9 


1.53 


1 


41.1 ±9.2 


> 0.03 



a 90% positional uncertainty radius, including the 0'.'64 systematic Chandra pointing uncertainty added in quadrature. 

b Binning in which the source was best detected (§ |2.2| 

c Background subtracted 

d Based on the fact that the hardness must lie in the range -1 to 1 

e See Tabled 
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Fig. 1.— ■ Hardness ratio as a function of number of 0.3 — lOkeV counts for the sources in 
Table [2] The four cases where there are possible matches to our VLA sources (Table [I]) are 
marked. 

Event filtering was performed using 3x3 pixel islands. This approach, compared to the 
use of 5 x 5 pixel islands, may result in a slightly larger fraction of cosmic rays being misiden- 
tified as real sources, but also has the benefit of being less likely to reject real astronomical 
sources. 

Hardness ratios were computed as HR = ((C 2 — Ci)/{C\ + C 2 )) where C\ is the counts 
(after background subtraction) in the 0.3 — 2keV band and C2 the counts in 2 — lOkeV. We 
tabulate HR in Table |H In some cases the formal errors on HR resulted in values outside 
the range — 1 < HR < 1; here we truncate the values at these limits and report upper or 
lower limits on HR as appropriate. In Fig. [TJ we plot HR as a function of number of counts 
for the sources in Table [21 

In Fig. [2j we show a multi-wavelength view of the field, created by combining the g-, 
R-, J-, H- and .£C.-band data from B07, and overlaying the B07 radio catalog (Table II]), and 
the X-ray catalog from Table [2] 



3. Multi-Frequency Matching 

3.1. Radio - X-Ray Catalog Matching 

We searched around the radio positions (Table [I]) for matches in the X-ray catalog 
(Table |2J). In order to determine the maximum radius where a match is believable, we ran 
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Fig. 2.— Optical / IR image of the field with radio and X-ray positions overlaid. The 
color image was made from the data from B07: Keck g- and i?-band images (in the blue and 
green channels respectively) and the mean of the PAIRITEL J-, H-, and i^ s -band images (in 
the red channel). Green circles are the radio sources from B07, as labelled in Table [IJ and 
white circles are the Chandra sources reported here, as labelled in Table [2j The sizes of the 
white ellipses correspond to the size of the Chandra sources as determined by the wavdetect 
software, except for cases where the size was 5" or smaller, in which case we show a circle 
with a radius of 5" for clarity. The green circles are plotted with a radius of 5" for clarity; 
typical radio position uncertainties are ~ 1" (Table [I]). The yellow dashed box shows the 
edge of the Chandra field. 
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Monte Carlo simulations. The areal density of Chandra sources in our image is ~ 1600 cleg -2 . 
We ran 1000 iterations of a routine where we created mock X-ray catalogs with the same 
density (over an area twice as large as our VLA field), and then searched for matches in 
this simulated catalog around the VLA positions. We calculated f(9), the fraction of VLA 
sources with a match in the catalog closer than radius 9, for both the mock and real Chandra 
catalogs. At small 9, we expect few matches in the mock catalog, but as 9 increases, the 
fraction of sources with matches also increases due to chance associations with neighboring 
sources. In Fig. pi we plot f(9) for both the mock and real catalogs. For the mock catalog, 
the match fraction is very closely approximated by the expected analytic form, 

fmock{9) = l-e-^ e \ (1) 

where p = 1.2 x 10 _4 arcsec~ 2 is the areal density of Chandra sources. For a match at a 
given radius 9, f m0 ck(9) gives an estimate of the probability that the match is due to chance. 
This assumes a uniform surface density of X-ray sources; if sources in the mock catalogs were 
clustered (as they are in the real catalog), f m ock{9) would tend to be lower over areas of the X- 
ray image away from concentrations of X-ray sources, leading to a tendency to overestimate 
the false match probability. Proper motion is also not accounted for; we discuss this further 
in Section H31 

For 9 = 11'.'7, f m0 ck(9) = 0.05, i.e., 5% of VLA sources have a match in the mock 
catalogs within this radius. If sources in the X-ray catalog are physically associated with 
sources in the radio catalog, we expect an excess of matches at small 9 over the expectation 
from the Monte Carlo simulations, and indeed this is seen when we compare f moc k{9) and 
freai{9) in Fig. [3j The 5% false match probability leads us to expect < 1 false match out 
of our 18 radio sources at a radius of 11''7, but for matches that are closer than this, the 
probability increases that matched sources are physically associated rather than superposed 
by chance. 

We used Equation [l] to calculate values of f m ock(9) (i. e., the probability that the associ- 
ation is due to chance) for radio sources with a match in the real Chandra catalog closer than 
ll'.' 7. Two of the steady sources, and two of the radio transients (one short, and one long) 



have possible Chandra counterparts, and we discuss each of these cases below (Sections 3.1.1 



and 3.1.2) 



Bandwidth smearing in the radio image is not a significant issue. At worst it increases 
the size of unresolved sources close to the edge of the field in the 5 GHz VLA A-array 
observations of B07 by around 70%, but in most cases is significantly less than this. Even in 
the worst case it should not make a large difference to the statistical significance of matches 
between the X-ray and radio catalogs. 
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Fig. 3.— Left: Fraction of sources, f(9), with a match at radius 9 m < 6, shown as a 
function of 9, where we matched the 18 sources in the B07 VLA catalog with our Chandra 
catalog (dashed line) and with 1000 Monte Carlo iterations of a mock catalog with the 
same areal density as our Chandra catalog (solid line). The curve from the Monte Carlo 
simulations is a very good match to the expected analytic form, f moc k(9) = 1 — e _7rpe , 
where p = 1.2 x 10~ 4 arcsec -2 is the areal density of Chandra sources. The real Chandra 
catalog shows an excess above random of matches to the VLA sources, showing that at least 
some of the radio sources are really associated with X-ray sources. Right: The same plot, 
but removing the 8 steady VLA sources and plotting the curves only for the 10 transient 
sources. An excess above random is still seen. 
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Due to the small number of matches between X-ray and radio sources, it is not possible 
to look for systematic offsets between the X-ray and radio astrometric frames. However, 
VLA astrometric accuracy is typically ~ Qf!l or better, so systematic errors are likely domi- 
nated by the 0'.'64 Chandra pointing uncertainty, which is already included in the positional 
uncertainties reported in Table [2j We show in § 3J3 that the Chandra data are well-matched 



to the optical data. Sub-arcsecond systematic offsets would have only a small effect on the 
confidence levels for the matches discussed below. 



3.1.1. Two X-Ray Matches to Steady Radio Sources 

Steady radio source 5S1 was matched to X-ray source X29 at a distance of lCf!79. 
The 90% positional uncertainty of X29 is 1'.'16. The positional uncertainties reported by 
B07 for source 5S1 evaluate to a 90% uncertainty of C/53. Adding these in quadrature 
gives a total 90% positional uncertainty of 1'.'28. Source X29 is extended, and 5S1 lies 
within the boundaries of the Chandra source (Fig. [2]). From Equation [IJ f(9 moc k = 10'.'79 ± 
1'.'28) = 0.0429^o;qq94, so such an association has a 4. 29^94% probability of happening by 
chance. There are 18 sources altogether (eight steady and ten transient), and the binomial 
probability of a false association for one or more of them is 54.6%, so the formal confidence 
for this association is just 45.4%. However, source X29 is bright and extended, and is clearly 
associated with a group of galaxies in the OIR image (Fig. [2]), one of which is clearly the 
host of the radio source, so the sources both do appear to have a physical association with 
the galaxy group (§ 3.3[ ). 



Steady source 5S5/8S2 was matched to source X58 at a distance of CfMO. The 90% 
positional uncertainty of X58 is 0'.'68. The 90% positional uncertainty of 5S5 is 0Y53. Adding 
these in quadrature gives a total 90% positional uncertainty of 0'. / 86, so the match is well 
within the uncertainties. The false match probability is 0.006to'oo6% 5 or 0-1% f° r a ^ least 
one out of 18 sources. So we can be 99.9% confident that the association is real. 

In Fig. [4] we show X-ray images with radio contours overlaid of the two fields with X-ray 
counterparts to our steady sources. 



3.1.2. Two Possible X-Ray Matches to Transient Radio Sources 

Long (duration ~ 2 months) transient source 8L1 was matched to source X124 at 
a distance of 1'.'70. The 90% position uncertainty for X124 is 0'.'95. The radio position 
uncertainties reported by B07 for source 8L1 evaluate to a 90% uncertainty of C'54. Adding 
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Fig. 4.— - Images of the two fields with X-ray counterparts to our steady radio sources 5S1 
(left) and 5S5 (right). The colors in the image represent the number of Chandra 0.3 — 10 keV 
counts per pixel. The image is binned in blocks of 4 x 4 pixels. The radio data are from 
the deep 5 GHz image, with contours from 15 - 35/zJy beam -1 in steps of 5/^Jy beam -1 for 
5S1, and from 50 - 200 /xJy beam -1 in steps of 50/xJy beam -1 for 5S5. The size of the blue 
ellipse for source X29 corresponds to the size of the Chandra source as determined by the 
wavdetect software. For source X58 the fitted size was < 5" and the circle is shown with a 
radius of 5" for clarity. 
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these in quadrature gives a total 90% positional uncertainty of 1'.'28. So the measured 
positional offset is just a little larger than the formal positional uncertainty. From Equation[TJ 
fmock{0 = 1"70 ± 1"28) = 0.001 l+aooio) so such an association has a 0.11±°;^% probability 
of happening by chance. There are ten radio transients, and the binomial probability of one 
or more of them being associated with a Chandra source by chance is 1.0%. Hence, we can 
be 99% confident that the association between 8L1 and X124 is real. We note, however, 
that X124 was detected with 8.8 ± 3.6 counts, so the X-ray source is not a very strong 
detection. Additionally, this source is not present in a catalog filtered using 5x5 pixel 



islands (Section 2.2). Examination of the event list before filtering shows that four out of 
21 events occurred at the same time in the same pixel. These events are likely not due to 
photons from the source, but may be due to a cosmic ray or hot pixel. The absence of X124 
from the 5x5 pixel filtered catalog means that we cannot rule out that it may be a spurious 
source due to cosmic rays or particle background. Additional X-ray data would be useful to 
confirm or reject this X-ray detection. 

Short duration transient source 5T7 was matched to source X19 at a distance of 5'.'57. 
The positional uncertainty of source X19 is C/83. The positional uncertainty of source 5T7 
is 3"06. Adding these in quadrature for a total uncertainty of 3'.'17, we see that the match 
radius is almost twice the formal uncertainties. However, the X-ray source is towards the edge 
of the X-ray field and appears somewhat extended, and the synthesized beam of the radio 
observations in which this transient was detected is relatively large (FWHM 28" x 13"); 
bandwidth smearing for this source is not a significant problem, though (broadening the 
radio source by around 10% in the radial direction). Equation [I] shows that the false match 
probability is 1.16^o'g5%. The binomial probability of one or more out of ten transients being 
associated with a Chandra source this close is 11.0%, so we can be 89.0% confident that this 
association is real. We discuss this source further in Section l3~5l 

In Fig. [5j we show X-ray images with radio contours overlaid of the two fields with 
possible X-ray counterparts to our radio transients. 



3.2. X-Ray Upper Limits at Radio Positions 

For the eight transients and three steady sources without a counterpart from our Chan- 
dra catalog within 11"7 (not including the three steady sources which are outside the area of 
our Chandra coverage), we measured counts from the Chandra 0.3 — lOkeV image. We used 
apertures with radii ranging from 5 - 10" (depending on how far away from the center of 
the Chandra image the positions were, to account for distortion of the PSF) centered at the 
radio positions. The counts are reported in Table [I] We also tabulate 90% Poisson upper 





40.00s 15h03m36.00s 



4.00s 14h59m40.00s 



Fig. 5.— Images of the two fields with possible X-ray counterparts to our radio transients 
8L1 (left) and 5T7 (right). The colors in the image represent the number of Chandra 0.3 — 
10 keV counts per pixel. The image is binned in blocks of 4 x 4 pixels. The radio contours are 
from 15 - 35/zJy beam -1 in steps of 5/zJy beam -1 for the two-month average VLA B-array 
data for 8L1, and from 150 - 250 /iJy beam -1 in steps of 50/iJy beam -1 for the single-epoch 
lower resolution VLA D-array data for 5T7. The size of the blue ellipse for source X19 
corresponds to the size of the Chandra source as determined by the wavdetect software. For 
source X124 the fitted size was < 5" and the circle is shown with a radius of 5" for clarity. 
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limits for the count rates following Gehrels (1986). 



3.3. Optical / IR matches 

We generated a catalog from the i?-band image of B07, who report an astrometric 
uncertainty relative to the International Celestial Reference System (ICRS) of 250 mas in 
each axis for this image. We matched this .R-band catalog to our Chandra catalog, to check for 
systematic astrometric shifts in our Chandra data. For the 40 Chandra sources with a match 
in the i?-band catalog within a radius of 10", we found a median offset of ARA = 0"27±3"91, 
ADec = 0'.'20 ± 1'.'60. We conclude that our Chandra data are astrometrically well-matched 
to the ICRS. We neglect these small offsets in our analysis of the matches, because they are 
smaller than the other errors. 

In Table [TJ we note the OIR matches as reported by B07 for the transient and long 
transient radio sources. We also report where there are OIR counterparts for the steady 
sources as seen in the LRIS and PAIRITEL images. Two sources (5S3 and 5S4) have clear 
counterparts which appear as faint, small angular size sources in the OIR images. One 
source (5S8) has a probable association with an IR source which is outside the footprint of 
the optical coverage. These three sources are likely AGNs or perhaps starburst galaxies at 
intermediate redshifts. None has an X-ray counterpart. 

Source 5S1 is associated with a galaxy in a group or cluster which appears to be at 
intermediate redshift. The cluster is coincident with an extended X-ray source (X29) which 



is presumably due to bremsstrahlung from the cluster halo (§ 4.8). Transient source 5T6 is 
also quite nearby, and may be associated with a galaxy at z = 0.25, as reported by BO 7, 
although as noted above the offset between the galaxy and the radio position is rather large 
to be plausible. It is possible that the potential host galaxy is a member (on the outskirts) 
of this cluster — it is 42" from the cluster radio source 5S1, and 35" from the centroid of the 
X-ray halo X29, which if the cluster is also at z — 0.25, correspond to 163 kpc and 136 kpc 
respectively, within the virial radius for even a poor cluster. 

For the sources with OIR identifications, we searched the NASA Extragalactic Databasdj 
for redshifts or other information about the counterparts, but no further information was 
found. 



http : //nedwww . ipac . caltech . edu/ 
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3.4. Proper Motions 

For nearby sources such as brown dwarfs, high proper motions may result in positional 
mismatches between the same source in the radio, OIR, and X-ray data, which can complicate 
the interpretation of the data. For the two radio transients with nearby X-ray matches, the 
time between the X-ray and radio observations was ~ 10 yr. From the offsets in Table [TJ 
this implies that if the X-ray and radio sources were exactly coincident at the time of the 
radio transient (and ignoring any systematic position errors), the progenitors of 8L1 and 5T7 
would have proper motions (with 90% confidence uncertainties) of /x = 0.18±0.12 arcsec yr -1 
and n = 0.50 ± 0.28 arcsec yr -1 , respectively. 



3.5. Source X19 

3.5.1. Spectral Fits 

Source X19 was detected with 40.5 counts in our X-ray data, and we are able to fit 
models to the X-ray spectrum (Fig. J6|, as opposed to simply predicting fluxes from basic 
models for the fainter sources. Using the XSPEC software, we fit the fluxes using three 
models: a power law, where the photon index and absorbing column were allowed to vary, 
and blackbody and bremsstrahlung models where the temperature and absorbing column 



were allowed to vary. Formal errors were determined using the Cash (1979) statistic. Best 
fit parameters with 90% confidence limits, along with the reduced x 2 values for the best- 
fitting models, and associated fluxes are shown in Table [3j 

The X19 spectrum is about equally well fit by an absorbed power law with T = 1.7 as 
by a bremsstrahlung model with kT = llkeV. A single-temperature blackbody is a worse 
fit to the data than these two models. There are insufficient data to perform two-component 

Table 3. Fits for Source X19 
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a 0.3 — lOkeV, unabsorbed, computed for the best fit parameters. 
b With three degrees of freedom. 
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Fig. 6.— Top: X-ray spectrum of source X19. The model shown is the best-fitting absorbed 
power-law with N H < 1.2 x 10 22 cm" 2 , and T = 1.7i{];jj (TableQ. Bottom: The Ax residuals 
(i.e. residuals divided by uncertainties), denoted x-> f° r the fit to the spectrum shown in the 
top panel. 
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fits to the spectrum. 



3.5.2. X-Ray Variability 

In Fig. [7j we show the X-ray light curve for XI 9. By eye there appears to be some 
variability but the error bars are large, and most of them are consistent with the weighted 
average of the Chandra counts. A Kolmogorov-Smirnov test shows a 32% chance that the 
source is constant, i.e., variability is detected at only the la level. 

The power spectrum for X19 is shown in Fig. [7} For most frequencies, the variability 
is consistent with the Poisson noise level with no additional variability from the source, 
although there is tentative evidence of variability at low frequencies. Additional epochs 
of X-ray data could help to determine if the source is truly variable on longer timescales. 
However, with only a small number of counts per second, it will remain difficult to detect 
variability in this source on short timescales using Chandra (unless it brightens significantly 
in the future). 



3.5.3. Multiwavelength Counterparts 

Although 5T7, the source apparently matched to X19, has no counterpart in the deep 
OIR imaging of B07 or O10, it is offset from the X19 position by 5'/ 57. On examination 
of the i^-band images of the field presented in O10 (their Fig. 1, bottom right panel), we 
noticed a faint source offset from the radio position, but which appears to be consistent with 
the X-ray position. The source is not present in the K s -band data of B07, so it must be 
fainter than the detection limit in B07 (K s ~ 18) but brighter than the detection limit in 
O10 (K a ~ 19.2). From its appearance in the O10 data, we estimate K s = 18.5 ± 0.3. The 
source is also not visible to H ~ 18.5 and J ~ 19.2 in the B07 data. 



We searched the Wide-field Infrared Survey Explorer (WISE) catalog (Wright et al. 



2010) for counterparts to X19, and found that it was detected in the two shortest wavelength 
bands (3.4 and 4.6 /zm, with magnitudes 17.16 ± 0.10 and 16.14 ± 0.13 respectively), and 
undetected at the two longer wavelengths (12 and 22 /im, to 12.9 and 9.1 mag respectively). 

On examination of the Second Palomar Observatory Sky Survey (POSS-II) digitized 
images ( Reid et aL]|1991 ) we also found faint sources at the same position in B- and i?-band, 



and no detection in /-band. Quoted limiting magnitudes are Bj = 22.5, R c = 20.8, and 
I c = 19.5. From the POSS-II images, we estimate Bj = 21.0 ± 0.5, R c = 20.0 ± 0.5, and 
I c > 19.5. 
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We also searched the Galaxy Evolution Explorer (GALEX; Martin et al. 2005) Release 



6 catalog for ultraviolet (UV) counterparts to X19, and discovered a source with a near- 
UV (1750 - 2800A) magnitude in the AB system of 22.43 ± 0.34, corresponding to a Vega 
magnitude of 20.66. In the far-UV (1350 - 1750 A), there is no GALEX detection to AB 
~ 22.3 mag, corresponding to a Vega magnitude fainter than 19.9. 

Postage stamp images at the position of X19 are shown in Fig. [8j 



3.5.4- Proper Motion 



As noted in Section 3.4, if the offset between the positions of 5T7 (seen at UT 1999 May 
4) and X19 (observed UT 2010 July 3) is due to movement of the source, a proper motion 
of 0.18 ± 0.12 arcsec yr™ 1 is required. However, the position of the OIR counterpart to X19 
as seen in the Bj, R c , and 7 c -band POSS-II data (taken at UT 1994 April 6, 1994 June 09, 
and 1997 June 28, respectively) is coincident with the WISE position (data taken in early 
2010), as well as the GALEX position (UT 2006 November 5). Therefore, proper motion 
cannot explain the 5'.'57 offset between the 5T7 and X19 positions, which are bracketed by 
epochs where the UV/OIR counterpart to X19 is not seen to move. 



3.5.5. Spectral Energy Distribution Fitting 

We performed spectral energy distribution (SED) fitting using the UV and OIR data 
for X19. A x 2 minimization using a library of 726 galaxy and QSO models, and 254 stellar 
models (including 100 brown dwarf models), was performed using the Le PHARE software 



flllbert et al.||2006[). The best-fi tting galaxy or QSO template, a QSOl from the SWIRE 
template library ( |Polletta et al.||2~0~0~7| redshifted to z = 1.29±{j;J5, had xl = 1-47; the best- 
fitting stellar model, a brown dwarf from Chabrier et al. (2000), was a poor match to the 



data (xl — 24). In Fig. [9j we show the SED of X19 along with the multiwavelength UV/OIR 
photometry. 

The X-ray light-curve and spectrum, and the positional coincidence of the X-ray and 
UV/OIR sources, make for a convincing argument that X19 is a QSO at z — 1.29. The 
offset between X19 and 5T7 must then be due to one of two possibilities: either 5T7 is not 
at z — 1.29, in which case X19 and its multiwavelength counterpart are unassociated with 
the radio transient, which remains unidentified in X-rays; or 5T7 is at z = 1.29, in which 
case it appears to be offset by 47 ± 27kpc from the X-ray position. We discuss this source 
further in Section Wjl 
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Fig. 7.— Left: Light curve for source X19 during our Chandra observation. The dashed 
line shows the weighted average of the Chandra counts. The dotted line shows the back- 



ground count rate. Right: Power spectrum for source X19, using the Leahy et al. (1983) 



normalization, where a power of 2 indicates consistency with Poisson noise with no additional 
variability from the source. 
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Fig. 8. — Postage stamp images of X19. From top left to bottom right, in order of increasing 
wavelength: Chandra; GALEX far-UV and near-UV; POSS-II Bj, R c , and I c \ J, H, and K s 
from B07; WISE channels 1-4. Overlaid on each image are radio contours are from 150 - 
250/iJy beam -1 in steps of 50/zJy beam -1 for the single-epoch VLA D-array data for 5T7. 
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Fig. 9.— Spectral Energy Distribution of X19. From left to right, in order of increasing 
wavelength, the data points are: GALEX far-UV and near-UV; POSS-II Bj, R c , and J c ; J, 
H, and K s from B07 and O10; WISE channels 1 - 4. The best-fitting QSO model, a QSOl 



from Polletta et al. (2007) at z = 1.29, is shown. 



26 



4. Transient Progenitors 

It is important to note here that it is unlikely that all the radio sources belong to the 
same population. Three, or maybe four of the steady sources appear to be associated with 
galaxies in the OIR images. These are likely AGNs or starburst galaxies. The remainder 
may be Galactic or extragalactic, but are likely distant or highly obscured. 

The transient sources with no counterparts may be very distant, highly obscured, or 
intrinsically faint in quiescence. Again, this is true whether they are Galactic or extragalactic, 
although in the latter case, as argued by O10, we would still expect to see more counterparts. 

The two transient sources with clear OIR counterparts, and perhaps the two with pos- 
sible counterparts, may, as noted by B07, be supernovae, although if so, they are unusually 
radio-luminous. For these four sources, no X-rays were seen in our Chandra image. This 
would not be particularly surprising if these sources are supernovae, particularly since the 
X-ray observations were taken between 11 and 27 years after the radio transients were seen. 



It has also been recently suggested (Nakar & Piran 2011) that the radio emission in the two 



cases with clear counterparts may be from neutron star - neutron star or neutron star - 
black hole mergers. 

Of the six transient sources with no OIR identification, two have possible X-ray coun- 
terparts in our Chandra image (although one is a marginal detection and the other is offset 
too much for the association to be certain), and four have only upper limits. The two sources 
with "possible" OIR IDs as noted by B07 (which we argue are also unidentified) also have 
only upper limits in our Chandra data. 

The count rates of the upper limits are similar to those for the detections, however, so 
it is conceivable that all the transients have similar count rates, of order 10~ 4 counts s _1 . 
They may or may not be due to a single class of progenitor. It is worth noting that the only 
short-duration transient with a possible X-ray counterpart (5T7/X19) is by far the brightest 
of the transients in the radio (Table Eh, which could indicate that it is a member of a different 
class of objects to the others. 

We can achieve an improvement in sensitivity by "stacking" the X-ray images at the 
radio positions of the six sources with no OIR identification and only X-ray upper limits 
(5T1, 5T3, 5T4, 5T5, 5T6, and 8T1; see Table [1|. Stacking is a well-established (e.g., 



Zibetti et al. 2005; de Vries et al. 2007) method for combining data from individual sources 
in order to study the statistical properties of objects which are undetected in single images. 
Prior information about positions at a frequency where sources are individually detected (in 
this case, the radio positions) is used to create "cut-out" images centered on those positions 
at a second frequency (here, from our X-ray data) where sources may not be individually 
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detected. If iV images each with rms a are combined, then the noise level of the stacked 
image will decrease approximately as a/y/N. This allows sources below the original detection 
threshold to be studied. Stacks can also be created at "blank sky" positions offset from those 
of the known sources, in order to determine the background noise level in the absence of any 
additional emission present at the source positions. 

We created a stack from the mean of six X-ray images aligned at the radio positions, and 
eight blank sky stacks; four at positions offset 30" in the cardinal directions, and four offset 
42" in the ordinal directions from the radio positions. In other words, we created stacked 
images at positions on a 3 x 3 grid, with spacing 30", centered on the radio positions. 

We measured Chandra counts in a 5"-radius aperture at the center of each stacked image. 
The stacks offset from the radio positions will be dominated by background photons, whereas 
the stack centered at the radio positions will contain any additional flux from the radio 
sources, if present (assuming proper motions smaller than a few arcseconds over the time 
between the VLA and Chandra observations). The mean and rms of the eight background 
values was 7.09 ± 0.89 counts. There were 7.93 counts in the aperture in the stack made at 
the radio positions. So there is less than la excess of X-ray counts at the radio positions of 
the non X-ray detected sources. 



4.1. Models, Fluxes, and Spectra 

We used the WebPIMMS v4.2 softwareSto convert the count rates into flux densities 
(for the two Chandra sources near to B07 transients) or upper limits (for the eight transients 
with no X-ray counterpart). The calculations are somewhat dependent on the model used 
for the X-ray emission. 



4-1.1. Simple Models 

We chose to model two cases: a neutron star, and a flare star. The neutron star case 
was modelled as a blackbody with energy kT = 0.1 keV, since isolated neutron stars older 



than ~ 10 y have temperatures < 0.1 keV (Yakovlev et al. 2005). We caution, however 



that the only isolated neutron stars which have been observed to date are relatively young, 
and may not be representative of older populations. Slow accretion onto magnetic neutron 
stars with fields of ~ 10 12 G can cause > 10% of the X-ray luminosity to be emitted from 



^http: //cxc .harvard. edu/toolkit/pimms . jsp 
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broadened cyclotron emission features in the hard X-ray (~ lOkeV) regime (Nelson et al. 



1993). Other authors, however (e.g. Zane et al. 2000) note that although the spectrum for 



such objects is slightly harder than thermal, it is not markedly different from a blackbody. 

For the flare star we used a bremsstrahlung model with kT = 1.0 keV. Both flare and 
post-flare spectra are often modelled in the literature using two-temperature models (e. g. 



Osten et al. 2010). Since we have no reason to believe that the progenitors were in outburst 



when observed with Chandra (in contrast to when they were seen to flare in our VLA 
observations), the post-flare model is appropriate here. We choose a 1 keV single temperature 



model, which is a valid approximation for both flare stars (Bildsten & Rutledge 2000) and 



brown dwarfs ( Imanishi et al.pOOl ; Preibisch et al.|[2005 ). Blackbody, Raymond-Smith, and 
thermal bremsstrahlung models give similar fluxes (to within 30%) in the 0.3 — lOkeV energy 
range. We choose a thermal bremsstrahlung model. 

More detailed spectral information could be used to constrain two-temperature fits or 
non-thermal components for both the neutron star and flare star cases, but since our sources 
are mostly undetected or barely detected, a simpler approach is warranted. 

In both the neutron star and flare star cases, we assumed a Galactic H I column density 
of 3 x 10 20 cm -2 (Kalberla et al. 2005). We show the resulting flux densities and limits in 
Table H 

For source X124/8L1, which was detected with only 8.8 counts, we have insufficient 
data to perform spectral modeling. However, we can use the HR constraints (0.38 ± 0.44) 
to constrain progenitor models. We used WebPIMMS to compute hardness ratios for black- 
bodies with a range of energies, and we find kT > 0.8 keV at 90% confidence. The best fit 
kT = 1.4 keV. 



4.2. Luminosities 

The X-ray luminosities of the transients (ignoring /c-corrections if they are extragalactic) 
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where So.3-iofceV is the model-dependent X-ray flux or upper limit from Table [4] or Table [3j 
and d is the luminosity distance to the progenitor in kpc. 



The radio luminosities of the transients are 
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where 5*5 qh z is the transient radio flux density at 5 GHz (Table El), or the upper limit (from 
B07) in quiescence (typically a few tens of /xJy), and d is the luminosity distance to the 
progenitor in kpc. 



4.3. Isolated Old Galactic Neutron Stars 



O10 suggest that the B07 transients may be due to isolated old Galactic neutron stars 
at distances between ~ 1 - 5 kpc. This implies X-ray luminosities around 10 29 - 10 31 erg s _1 
for X124, and similar values as upper limits for the X-ray non-detections. This is consistent 
with the typical X-ray luminosities, < 10 31 erg s _1 , of old neutron stars accreting from the 
interstellar medium (ISM; Treves et al.|2000 Ostriker et al.fl~9~7"0 ). This luminosity depends 
on the accretion rate from the ISM, which in turn depends on the ISM density and the 
neutron star space velocity, and also on the field strength and spin rate of the neutron star. 
It also depends on the details of Bondi accretion in the presence of magnetic fields, which 



are not well understood (Toropina et al. 2003 2005 Arons h Lea 1976 1980). 



However, although the measured X-ray fluxes and upper limits are consistent with a 
population of isolated neutron stars at distances between ~ 1 - 5 kpc, the hardness ratio 



measurements for X124 (Section 4.1.1) exclude a simple one-component model with kT 



0.1 keV, the expected energy for neutron stars. Since X19 is clearly not a neutron star, 
and X124 is the only other source for which we can compute hardness ratios, we cannot 
extrapolate this conclusion to the remainder of the X-ray sources. Additionally, as previously 
noted, accreting neutron stars may also have non-thermal components, and the spectra of 
old isolated neutron stars may not be identical to that of younger cooling neutron stars. 
Also, as previously noted, X124 may in fact be a spurious detection. 



Rutledge et al. (2003) note that optical luminosities of isolated neutron stars are around 

5 orders of magnitude fainter in the optical than in the X-ray. Our X-ray fluxes of < 10~ 14 erg 

s _1 cm~ 2 would correspond to < 10~ 19 erg s _1 cm~ 2 optical fluxes for this model, or R ~ 35 

- out of reach even for the planned James Webb Space Telescope. Accretion can cause 



brightening by a factor of a few in the optical (Zane et al. 2000), but not likely enough to 
make our objects detectable. 



4.4. Magnetars 



Another class of potential neutron star progenitors discussed by B07 and O10 are soft 
gamma-ray repeaters (SGRs) or magnetars. These are attractive as potential progenitors 
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because they are known to exhibit dramatic flares (Gaensler et al. 2005), and would have 



no OIR counterparts in quiescence in our data. The X-ray emission from magnetars almost 
always has a blackbody component with kT ~ 0.5 keV, and sometimes also has a power-law 



component (Mereghetti 2008 ). Our X-ray data for X124 are not sensitive enough to constrain 



multi-component fits. 

Typical magnetar X-ray luminosities are ~ 10 33 erg s _1 , so our sources would have to be 
at > 30 kpc, or underluminous if closer, to explain their non-detection or marginal detection 
in X-rays. As noted by B07 and O10, the surface density of magnetars is too low to explain 
all the B07 transients as being due to this class of object. 



4.5. Flare Stars 

Low mass stars and brown dwarfs are numerous, optically faint, and active at radio 
wavelengths in 5 to 30% of cases ( Bergerp006 ) . M-dwarfs at distances as large as lkpc have 



been identified as the dominant contribution to optical transient rates (Becker et al. 2004 



Kulkarnifc Rau 2006). 



In an extreme case, the dMe star EV Lac was observed to undergo a factor of ~ 150 
flare to a radio luminosity ~ 10 15 erg s _1 Hz" 1 , with a rise time of minutes and a decay time 



of hours (Osten et al. 2005). Flare star (i.e., dMe star) radio luminosities in quiescence are 

10 15 erg s _1 Hz _1 (Giidel 2002), implying even extreme flares 
10 17 ' 



generally no brighter than 
ought to be no brighter than 



erg s Hz - . The radio luminosities of our transients 
(Equation |3j Table [2| are L 5Giiz ~ 10 18 (<i/l kpc) 2 erg s _1 Hz _1 . So if the progenitors are 
flare stars, they must be closer than 1 kpc, or else the transients must be due to unusually 
dramatic events. 



From the X-ray - radio relation (Giidel 2002), our quiescent X-ray flux limits imply 



quiescent radio flux densities of < 1 /iJy (consistent with the non-detection of the transients 
in the deep radio image of B07 to depths of a few tens of /i Jy) . The inferred X-ray luminosities 
(Equation [2j are < 10 30 erg s _1 for progenitors at lkpc, plausible for flare stars (Giidel 



2002). Predicting X-ray luminosities during outburst from the peak radio luminosities of 



the transients (assuming the X-ray - radio relation holds during outbursts), we obtain X-ray 
luminosities of Lx ~ 10 32 (cf/l kpc) 2 erg s . 

However, the non-detection of the transients to K = 20.4 mag by O10 implies that if 



they are due to M5 stars, they must be at d > 1200 pc (Patten et al. 2006). So unless the 



radio flares are unusually luminous, this implies that if flare stars are the progenitors of the 
transients, they must be of later spectral type than M5. 
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010 suggest that brown dwarfs at typical distances of around 200 pc are plausible pro- 
genitors for the B07 transients, although they note that we would expect the radio emission 
to be strongly circularly polarized, in contrast to the low levels of polarization observed 
by BO 7. If they have spectral types as late as T8, they would be too faint to be detected 



by O10 if at distances > 50 pc (Patten et al. 2006). Our X-ray fluxes imply luminosities of 



< 10 27 erg s x at 30 pc. This is a typical luminosity for quiescent emission from brown dwarfs 



(Preibisch et al. 2005); sources with non-detections in our X-ray data could be intrinsically 



less luminous, more obscured, or further away. 



4.6. Proper Motion of 8L1 (RT 20010331) / X124 



As noted in Section [3.5.4[ the offset between X19 and 5T7 cannot be explained by proper 
motion. There is no UV/OIR counterpart to X124, however, so it is plausible that a proper 
motion of /i = 0.18 ± 0.12 arcsec yr _1 is responsible for the 1'.'70 ± l'/09 offset between X124 
and 8L1. 

We can use the 90% confidence limits for the proper motion, along with plausible 
ranges of tangential velocity, v t = 4.74/irfkm s"" 1 , for different classes of progenitors, to 
infer plausible ranges of distance, d (in parsecs) to X124/8L1. A typical low mass star with 
Vt < 20 km s~ would have to be closer than ~ 60 pc. Some dwarf stars have tangential 



velocities > 100km s , however (Faherty et al. 2009), so X124/8L1 could plausibly be a 



high proper motion object at < 300 pc. This would appear to marginally rule out an M star 
progenitor for this source (since, as discussed above, M stars should have been seen by O10 
unless at > lkpc), but the proper motion is consistent with a brown dwarf progenitor. 

The measured proper motion is also not unreasonable for a neutron star at ~ 100 pc 
(Ofek 2009), particularly if the source is in the high- velocity tail of the distribution. A 



neutron star with tangential velocity v t < 3000 km s 1 (Cordes & Chernoff 1998 Faucher- 



Giguere & Kaspi 2006) could explain the proper motion of X124/8L1 if at < 9kpc. The 
lower limit for distances of old isolated neutron star progenitors from the OIR non-detections 
(> lkpc), as noted above, implies that the tangential velocity must be > 330km s~ . 

For both the flare star and neutron star cases, relatively small systematic position offsets 
could reduce the size of the measured proper motion, bringing upper limits on the distance 
from transverse velocities into better agreement with lower limits from to the K-b&nd non- 
detections by O10. If the faint X-ray detection is not in fact due to an astronomical source, 
we cannot constrain the proper motion of the progenitor. 
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4.7. Progenitor of 5T7 (RT 19990504) 



The X-ray luminosity of X19 is 2 x 10 44 erg s 1 , typical for X-ray selected AGNs (Lusso 



et al. 2010) 



The non-detection of X19 to S$ghz = 0.117 mJy in the deep radio image of B07 
implies that the quiescent radio luminosity of the quasar is L 5 ghz < 10 31 erg s _1 Hz -1 . 
Radio loud AGNs are defined as those having radio to optical luminosity ratio R > 10, 



where R = L 5 ghz/Lb (Kellermann et al. 1989). Our estimate of the B-band flux density 



Sb ~ 0.016 mJy, implies R < 7 from the deep radio data, so the AGN is not radio loud in 
quiescence. 

As noted above, the facts that X19 appears somewhat extended due to off-axis broad- 
ening of the Chandra PSF, and that the VLA D-array radio beam is rather large for 5T7, 
suggest that the uncertainty in the position may be somewhat higher than the formal un- 
certainties. If, in fact, 5T7 is associated with the AGN, this would imply that the AGN 
brightened by a factor of > 10 (see B07, Figure 4) to L 5 ghz ~ 7 x 10 32 erg s _1 Hz _1 , and 
then faded back by a factor of > 10 into undetectability on a timescale of days (and been 
brighter when detected as a transient by a factor > 60 relative to its mean luminosity on a 
timescale of 22 years; see B07, Table 2). During the flare, R ~ 440, so the AGN would have 
changed from radio-quiet to radio-loud in this scenario. 

It is more likely, however, due to the positional offset, that the radio source is unassoci- 
ated with the nuclear X-ray source. If the radio transient is also at z = 1.29, the positional 
offset implies it is 47±27kpc from the nucleus. In this case, the L 5 ghz ~ 7x 10 32 erg s _1 Hz -1 
radio transient could plausibly be due to a GRB afterglow in the outer regions or halo of the 



AGN host galaxy, although it would be among the brightest yet seen (Frail 2011). 



If in fact 5T7 is a chance association with X19, then the radio transient remains unde- 
tected in X-rays and should be considered along with the other non-detections as a possible 
brown dwarf or neutron star. 



4.8. Galaxy Cluster 

For source 5S1, apparently associated with a galaxy cluster, we used WebPIMMS to 
model the X-ray source X29 as a power law with T = 1 (optically thin bremsstrahlung) at 
z = 0.25, and a Galactic H I column density of 3 x 10 20 cm -2 , which results in a luminosity 



(extrapolated to rest frame 0.05 - 40keV) of ~ 5 x 10 43 ergs 1 . Following Reiprich & 
Bohringer (2002), we infer a mass of ~ 3 x 10 14 M for the cluster. The cluster could be part 
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of a structure associated with Abell2047 (see B07), which is centered 10'6 away, or ~ 2Mpc 
if both are at z ~ 0.25, perhaps explaining the fact that several of the B07 sources seem to 
be at this redshift. 



5. Conclusions 

Of the ten B07 transients, eight are undetected in X-rays; one is associated with an 
X-ray source of marginal (2.4a) significance, which may be due to cosmic rays or particle 
background; and one may be associated with a QSO or its host galaxy at z — 1.29, or 
may be a spurious positional coincidence. Because of the uncertainty regarding these two 
associations, it is plausible that all of the B07 transients remain undetected in our X-ray 
images. 

The measured X-ray flux for the transient with a marginal X-ray detection (X124/8L1), 
and the upper limits for the undetected transients, are consistent with the X-ray luminosity 
expected for brown dwarfs at distances ~ 100 pc, one of the possible progenitors discussed 
by O10. The small offset between the X124 and 8L1 positions could plausibly be due to the 
proper motion of such an object. 

If instead the transients are due to flaring M stars, most of the radio transients ought to 
have been detected in quiescence in both the X-ray and OIR images, unless the flares were 
extreme events (not out of the question, since we know the radio flux densities increased by 
a factor > 100) from progenitors at distances > lkpc. 

The X-ray fluxes are also consistent with the suggestion of O10 that the transients 
could be isolated neutron stars. The radio and OIR fluxes are consistent with this picture 
as well. If X124 is truly a counterpart to 8L1, hardness ratio constraints appear to rule out 
a simple one-component spectral model, but there are insufficient data to attempt a more 
sophisticated fit with additional thermal or non-thermal components. 

It is possible that the B07 transients are a mixture of some of the above classes of object, 
or some of them may be due to some unknown class of transient. 

If the progenitors have X-ray fluxes slightly fainter than our sensitivity limit, deeper 
observations would be expected to detect more of the transients in X-rays, as well as im- 
proving the significance of the fainter detection (or rejecting it as a real source). X-ray or 
radio monitoring campaigns lasting for several months would have a good chance of catching 
these sources in outburst since for the most likely classes of progenitor, flares are likely to 
repeat. Near- or mid-infrared data would also be useful to confirm or refute the brown dwarf 
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hypothesis. 
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Table 4. X-ray Flux Densities for Radio Transients 



VLA ID 


Chandra 


0.3- lOkeV rate 


5o.3-10feeV 


(NS) a 


So.3- 


-lOkeV 


(FS) b 


$5 GHz C 




ID 


(counts s _1 ) 


(10~ 15 erg s~ 


1 cm~ 2 ) 


(io- 15 


erg s — ] 


cm -2 ) 


(mJy) 


5T1 (RT 19840502) 




< 2.00 x 10" 4 


< 6.17 




< 1.92 






0.45 


5T2 (RT 19840613) 




< 2.51 x 10~ 4 


< 7.74 




< 2.41 






0.57 


5T3 (RT 19860115) 




< 1.09 x 10~ 4 


< 3.36 




< 1.05 






0.37 


5T4 (RT 19860122) 




< 2.00 x 10~ 4 


< 6.17 




< 1.92 






1.59 


5T5 (RT 19920826) 




< 1.17 x 10~ 4 


< 3.61 




< 1.12 






0.64 


5T6 (RT 19970528) 




< 2.67 x 10~ 4 


< 8.23 




< 2.57 






1.73 


5T7 (RT 19990504) 


X19 


8.20 x 10~ 4 


25.3 d 




7.88 d 






7.04 


8T1 (RT 19970205) 




< 1.76 x 10~ 4 


< 5.43 




< 1.69 






2.23 


5L1 (RT 19870422) 




< 4.29 x 10~ 4 


< 13.2 




<4.12 






0.51 


8L1 (RT 20010331) 


X124 


1.78 x 10" 4 


5.49 




1.71 






0.70 



a Unobscured flux for a basic neutron star model (blackbody with kT = 0.1 keV) 

b Unobscured flux for a basic flare star / brown dwarf model (thermal bremsstrahlung with kT = 1.0 keV) 

c From column 4 of Tables 2 and 3 of B07. For the 8 GHz transients, we assume a flat spectral index between 5 and 
8 GHz, consistent with the discussion in O10. 

d The X-ray source appears to be an AGN, and values from more sophisticated spectral fits for this source from Table|3J 
(see also Section |3, 5. 1\ should be used. The X-ray source may not be physically associated with 5T7 (Section|3.5.5[l. 



